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Abstract. Rotation curves are often used to estimate the mass distribution of spiral galaxies, assuming that the 
circular velocities of the interstellar medium balance with the galactic centrifugal force. However, non-circular 
motions caused by a non-axisymmetric gravitational potential, such as a stellar bar, may disturb the velocity 
field, resulting in errors in mass estimation, especially in the central regions of galaxies. This is because the line- 
of-sight velocity depends on the viewing angles in a non-axisymmetric flow. Observing rotation curves of edge-on 
galaxies in time-dependent numerical simulations from different viewing angles, we obtain errors in the estimation 
of galactic mass from the rotation curves. In the most extreme case, the ellipticity of gas orbits is as high as e ~ 0.8 
in the central regions, even if the bar potential is weak. When rotation curves are defined as the highest velocity 
envelope of position-velocity diagrams, the mass estimated from the rotation curves is larger than the true mass 
by a factor of five for 15% of the viewing angles, and the ratio between the apparent mass and true mass is less 
than six for any viewing angle. The overestimation in mass occurs more frequently than the underestimation. 
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1. Introduction 

Rotation curves are a major tool for determining mass 
distribution in spiral galaxies. Assuming spherical mass 
distribution, the galactic mass M within a galactocentric 
radius R is estimated by 



M = RV 2 /G, 



(1) 



where V and G stand for a circular velocity and gravita- 
tional constant. 

Rotation curves in their outer regions are generally 



about 10 Mq within a central 100 pc radius (Sofue, 1996 
|Sofue et al.| , |1999| ; [Takamiya fc Sofuc| , |2000| ). 



However, the gas in galactic disks does not necessarily 
show pure circular rotation, especially in the central re- 
gions. Bar-like distortions of the stellar system can drive 
non-circular (elliptical) motions for the gas. As a result, 
the apparent rotation curves do not represent the correct 
mass distribution. For example, if the elliptical orbits were 
aligned by chance with the linc-of-sight, we would overes- 



timate the mass (Sakamoto et al., 1999). The effect of non 



circular motion on position-velocity diagrams has been 



flat ( 


Rubin et al., 


198C, 1982, 


1985), indicating massive intensively studied in theoretical calculations 


^Bureau & 


dark hilos surrounding their optical disks (Kent, 


1987). Athanassoula , 


1999]; Athanassoula & Burcar: 


, ] 


1999 


)• 



In their central regions, many rotation curves rise steeply 
from the centers, reaching the high velocity seen in the 
outer flat rotation curves: typically 100—300 km s _1 within 
a central 100 pc radius ( Some , 1996| ; Sofue et al. , 1999| ). 



These high velocities may indicate central massive cores of 



In this paper, we quantitatively study the errors in es- 
timating the mass from rotation curves in galaxies with 
a weak bar, and calculate the probability that the ob- 
served mass suffers from such errors. Even if the bar- 
like distortion of the gravitational potential is very weak, 
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the gas velocity-field can be non-axisymmetric (Wada, 
1994). In order to obtain the velocity-field of the gas in 
a weak bar potential, we performed Smoothed Particle 
Hydrodynamics (SPH) simulations (Sect. ||). Using the 
numerical results, we estimate the probability of overesti- 
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mating the central galactic mass (Sect. ||). We discuss the 
implications of our results in Sect. [|. 

2. Numerical Calculations 

2.1. 2-D SPH Methods 



We perform two-dimensional SPH calculations (Lucy 



1977; Gingold & Monaghar, 1977), and reproduce gas 



motions in a weak stellar bar potential. We adopt the 



SPH formulation by Bcnz (1990), and use the spline kernel 



(Monaghan & Lattanzio, 1985) with the modification for 



its gradient ( Thomas fc Couchman , 1992). The correction 
term for viscosity flBalsara , 1995|) is taken into account 
to avoid large entropy generation in pure shear flows. The 
SPH smoothing length h varies in space and time, keeping 
the number of particles within the radius 2h at an almost 
constant of 32 according to the method of Hcrnquist fc 
Katz ( 1989 ). The leapfrog integrator is adopted to update 
positions and velocities. We use 3 x 10 4 particles to rep- 
resent the gas disk. 

2.2. Galaxy Disk and Weak Bar Potentials 

We take the g alax y potential for a wea k bar used in 



Sanders] (|1977|) and [Wada & Habe| (|1992| ); the oval po- 
tential for a barred galaxy $>(R,8) is represented in ax- 
isymmetric and bisymmetric parts by 



$(i?, 9) = <!> (R) + $ 6 (i?) cos 26», 



(2) 



where the first term is an axisymmetric potential for a 
disk, 



$ (i?) 



-C 



1 



(ff+a 2 ) 1 / 2 ' 



(3) 



and the second term is a bisymmetric one for a weak bar 
aR 2 



$ fc = -e C 



(R 2 + a 2 ) 2 ' 



(4) 



■«"max IS Cal- 



The normalization coefficient C = (27/4) 1 / 2 c 
culated from the core radius a and maximum rotation ve- 
locity w max - £o is a parameter to control the strength of the 
bar. R and 9 stand for the galactocentric radius and az- 
imuthal angle from the bar respectively. We are interested 
in the notion that noncircular motions in a bar appar- 
ently indicate a central massive component. We thus do 
not consider any massive central component in our poten- 
tial, such as a bulge ( Athanassoula fc Burcau| , |1999| ) or a 



massive black hole (Fukuda et al., 1998, 2000) 



Our potential model has the benefit of being capable 
of analytically investigating gaseous orbits in the bar po- 
tential (Wada, 1994), and has been well-studied in numer- 
ical simulations for the bar-driven gas fueling into galactic 
centers (Wada & Habe, 1992, 1995), the gas kinematics in 



the Galaxy (Wada et al., 1994), the spatial distribution 



of mass-to-light ratio in a galaxy NGC 4321 (Wada et al 



199S), and the effects of a central black hole (Fukuda et 
al.[ [1998| , p000| ). 




4 6 
R [kpc] 



Fig. 1. Rotation curve from the axisymmetric potential, 
eq. (§), with the core radius a — \/2kpc and maximum 
rotation velocity i> max = 220 km s -1 . 




R [kpc] 



Fig. 2. Radial changes of frequencies, Sl(i?) and il(R) ± 
k(R)/2. Four horizontal lines represent the pattern speeds 
of the bar, i.e. 0.4, 0.8, 1.0, and 1.5 X (Q - K/2) max . 



We fix a — \/2kpc and v max = 220 km s -1 ; the corre- 
sponding rotation curve is shown in Fig. |l|. The gas reaches 
the maximum circular rotation velocity at R = 2 kpc 
with a rotational period of 56 Myr. Fig. || shows the ra- 
dial changes of frequencies, SI and 51 ± k/2, where and 
k are circular and epicyclic frequencies respectively. We 
set the pattern speed of the bar fib at 0.4, 0.8, and 1.5 
times the maximum of f2 — k/2, indicated by horizontal 
lines. Models with Ob = 0.4 and 0.8 x (O — /-t/2) max have 
two inner Lindblad resonances (ILR), while those with 
fib = 1-5 x (51 — /i/2) max have no ILR. eo is set to 0.05, 
0.10, and 0.15. Our nine models are listed in Table [j]. 

2.3. Initial Conditions 

The gas is initially distributed in a uniform-density disk 
with an 8 kpc radius, following pure circular-rotation that 
balances the centrifugal force. The gas temperature is as- 
sumed to be a constant 10 4 K, corresponding to the sound 
speed of about lOkrns -1 , throughout evolution. The to- 
tal gas mass is assumed to be 5% of the total stellar mass 
within the radius of 8 kpc. The results are not significantly 
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model fib eo ILRs CR OLR (/9b/po)max T bar 

(kms-'kpe" 1 ) (fi- K /2) max (kpc) (kpc) (kpc) (kpc) (%) (Myr) 

A 49 1.5 0.05 -- 4.0 5.4 5 130 

B 49 1.5 0.10 -- 4.0 5.4 10 130 

C 49 1.5 0.15 -- 4.0 5.4 16 130 

D 26 0.8 0.05 1.1 3.1 6.2 8.3 5 240 

E 26 0.8 0.10 1.1 3.1 6.2 8.3 10 240 

F. .... 26 0.8 0.15 1.1 3.1 6.2 8.3 16 240 

G 13 0.4 0.05 0.6 5.9 10.0 13.3 5 470 

H 13 0.4 0.10 0.6 5.9 10.0 13.3 10 470 

I 13 0.4 0.15 0.6 5.9 10.0 13.3 16 470 



Two free parameters, i.e. bar pattern speed fib and bar strength eo, and radii of inner Lindblad resonances (ILR), corotation 
resonance (CR) and outer Lindblad resonance (OLR), maximum density ratio of the bar over the disk (pi,/po)max, and 
rotational time-scale of the bar Tb ar [= 27r/fib] at R = 2 kpc are tabulated. 



affected by the total gas mass, because thermal pressure 
is much smaller than the rotational energy, and we do not 
calculate self- gravity of the gas. We advance the calcula- 
tions up to about 500 Myr. 

2.4. Gas Dynamical Evolution 

Gas dynamics in a barred potential have been well studied 



in numerical simulations I 


Wada & Habe 




1992; 


Heller & 


Shlosman 




1994; 


Pincr et al. 




1999; 


Fukuda ct al. 




1998; 


Athanassoula & Bureau 




1999 


). Our models evolve 


con- 



sistently with these simulations. Fig. ||, model E, shows 
a typical evolution. Three phases of the evolution can 
be seen in this model: (a) linear perturbation phase, 
t ~ - 50 Myr, (b) transient phase, t ~ 50 - 250 Myr, 
and (c) quasi-steady phase, t > 250 Myr. 

The characteristic structure appearing during the evo- 
lution depends strongly on the positions of resonances, 
i.e. the pattern speed of the bar f^. In phase (a), lead- 
ing and trailing spiral arms are formed around the inner 
(R = 1.1 kpc) and outer (3.1 kpc) ILRs respectively at 
t — 36 Myr. Thes e resonant-dr iven spirals are expected 
in a linear theory (Wada, 1994| ). While the outer trailing 
arms remain with increasing density contrasts, the inner 
leading arms evolve into an oval ring, or a gaseous bar 
(t = 71 — 107 Myr), i.e. phase (b). The oval ring first leads 
the stellar bar (71 Myr), rotating opposite to the gas ro- 
tation (107 Myr), and being aligned with the stellar bar 
(250 Myr), and thereafter, the system develops toward a 
quasi-steady phase, i.e. phase (c). The ellipticity of the 
nuclear ring grows as high as e ~ 0.8. The ripple seen 
in the outer arms at t = 394 Myr would originate in the 

Gas dy- 



1995) 



Kelvin-Helmholtz instability (Piner ct al. 
namics and structure in the inner region of the disk are 
not affected by this instability. 

Fig. U and || display the final snapshots and velocity 
fields for nine models. Different pattern speeds fib and 
bar strength eo are arranged vertically and horizontally, 



respectively. It is evident that the final structure depends 
strongly on f2b, while eo changes only the density con- 
trasts. Model A, B, and C have no ILRs, thus no spiral 
arms or ring in their inner regions are formed. The outer 
spiral arms are formed outside the radius of the corota- 
tion resonance (CR) due to the outer Lindblad resonance 
(OLR). Model D and F resemble model E. Model G, H 
and I also have arms and rings similar to those in model 
E, but at different radii, corresponding to the location of 
the ILRs. Fig. |5| clearly show that, in models D, E, F and 
G, most gaseous orbits are aii-like, while in model H and 
I, the large separation and low density between the two 
ILRs suffice to leave the gases on X2-orbits, which form a 
stable oval ring, nearly perpendicular to the stellar bar. 



3. Galactic Mass Derived from Rotation Curves 

In this section, we compare the apparent rotation curves 



obtained from the numerical results (Sect. 3.1) and the 
true rotation curves, then we get probability to overes- 
timate/underestimate the galactic mass (Sects. 3.2 and 
3.3). 

3.1. Observing Rotation Curves in Models 



We obtain a position-velocity (p-v) diagram by observing 
our calculated gas disks edge-on, and then determine a ro- 
tation curve from the p-v diagram. We assign gas particles 
in a position-velocity grid using the cloud-in-cell method 
(Hockney & Eastwood, 1981); the spacing for the grid is 
set to a typical resolution in recent interferometry obser- 
vations of the CO gas for Virgo galaxies, i.e. 100 pc (~ 1") 
in space and 5 kms" 1 in velocity (Sofue in private commu- 
nication). Then we determine a rotation curve by tracing 
the gas at the highest velocity for each radius in the p-v 
diagram. Some examples for p-v diagrams and rotation 
curves are shown in Fig. ^ 
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Fig. 3. Gas dynamical evolution in model E: (fib, £o) = (0-8, 0.10). The stellar bar runs horizontally, and the gas 
rotates counterclockwise. After t ~ 250 Myr the system reaches a state of quasi-equilibrium. 



There are several ways to derive a rotation curve from 
an observed p-v diagram. One traces the peak-intensity 
velocity or intensity-weighted mean velocity at each ra- 



dius ( 


Rubin et al., 


1980, 


1982 


, 1985; 


Mathewson et al. 


1992; 


Mathewson & Ford, 


i99e 


) , while another traces the 



20% envelope of the peak-intensity velocity at each radius 
( jBofut , 1996). These intensity-based methods cannot be 
applied to our density-based p-v diagram, because inten- 
sity is not a simple function of density, especially in edge- 



on systems. Compared with these methods, our method 
provides generally higher velocity. 



3.2. Error Estimation in Rotation Curves and Mass 

We compare the observed rotation curves V bs(R) derived 
from p-v diagrams in simulations (Sect. 3.1), with the 
true rotation curves V^,ot(-R) from the gravitational po- 
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Fig. 4. Final snapshots for all nine models. Different pattern speeds fib and bar strengths eo are arranged vertically 
and horizontally, respectively. The stellar bar runs horizontally, and the gas rotates counterclockwise. 



tential (&o(R)). We estimate the errors in rotation curves 
by defining a function, i.e. 

V oha (R) 



ot(R) 



V pot (R) : 



(5) 



and examining whether a exceeds an arbitrary critical 
value a cr it as 



a(R) > Ocrif 



(G) 



The ratio of the mass M b s , derived from an observed 
rotation curve V b s , over the true mass M pot from the 
potential is, 



T (iJ)E^[=a(fl) 2 ]; 



pot 



(7) 



the second equal comes from Eq. 
comes equivalent to 

-y(R) > 7crit, 



M. Then Eq. 



be- 



(8) 



crit ' 



where 7 cri t = oc, 

A rotation curve rises steeply from a galactic cen- 
ter, having a peak, or at least a shoulder, at an in- 
nermost region, then reaching the flat rotation. A cen- 
tral galactic mass is always estimated at the radius of 



the peak or shoulder in observations (see Sofue et al 



1999). We thus consider the case that the mass is over- 



estimated/underestimated at the radius of the first peak 
or shoulder. The radius depends on, and changes with a 
viewing angle for the gas disk (see Fig. ph. Hence, we de- 
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0/05 0.10 0J5 e 

Fig. 5. Final velocity fields for all nine models. The arrangement is the same as in Fig. ||. Arrows are drawn for 1500 
out of total 30000 gas particles. 



fine a reference region of 400 pc < R < 2 kpc [] to which 
the above criterion , i.e. Eq. (|^), is applied; in our models 
the first peak or shoulder always fall in this region. If Eq. 
@ or (j|) is satisfied in this reference region, our observed 
rotation velocity differs from the true velocity at least by 
a factor of a C rit- 

3.3. Probability of Erroneous Mass Estimation 

We define the probability that the central rotation veloc- 
ity is larger/smaller by a factor of a cr it, as the fraction 
of viewing angles. We obtain rotation curves by observing 



1 The absolute values are not important because our results 
are almost scalable for an arbitrary core radius. 



the gas disk for viewing angles at ten-degree intervals, and 
calculate the probabilities for three-hundred snapshots in 
each simulation. We hereafter describe the probability as 
P[a > ttcrit] or Pa(ce C rit)- Fig. |?| shows a time evolution 
of P a in model E. Corresponding to the three evolution- 
ary phases (Sect. 2.4), the probability rises steeply in the 
linear perturbation phase, approaching a constant value 
in the transient phase, then reaching the constant in the 
quasi-steady phase. We also calculate the averaged prob- 
ability in the quasi-steady phase, i.e. 300 — 500 Myr, and 
describe it as P a . Similarly P 7 (7 cr it), the probability that 
the estimated mass differs from the true galactic mass by a 
factor of 7crit, is defined. P 1 is a time-average of P 7 (7 cr it)- 
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Fig. 6. Position-velocity diagrams and rotation curves in model E, observed from a variety of viewing angles, at 
t = 500 Myr. Solid lines show the distorted rotation curves derived from the calculated p-v diagrams (gray scale), 
while dashed lines indicate the rotation curve derived from the model potential. The mark "Detected" indicates that 
our criterion with a cr it = 2 is satisfied. 




500 



Myr 



Fig. 7. Change of the probability P[a > a CI n] = P a (a C iit) 
that Eq. H will be satisfied. a cr it = 1.6 and 2.0 indicate 
the cases that an observed rotation velocity exceeds the 
true value by a factor ofl.6 and 2.0 times, respectively. 



Fig. H shows the averaged probability P a as a func- 
tion of a cr ;t for all nine models. Different pattern speeds 
fib and bar strengths eo are arranged vertically and hori- 
zontally, respectively, as in Fig. || If the gas follows pure 
circular rotation, i.e. V Q h s = V po t, then these P a -profiles 
must be a step function: P a — 1 for a cr ;t < 1 and P a = 
for cc cr it > 1. However, the non-circular motion changes 
the P a -profiles as seen in the plots. In model E for ex- 
ample, P a — 0.8 at a cr it = 1.0 means that 80% of the 
rotation curves Vobs observed from random angles would 
apparently show higher velocities than the true rotation 
curve Vp t which traces the mass, and another 20% would 
show lower velocities than the true one. Thus the overes- 
timation in mass occurs more frequently than the under- 
estimation in model E. P a — 0.0 at a cr it = 2.4 indicates 
that the observed rotation curve cannot be overestimated 
by more than the factor of 2.4 in model E. 

We discussed in Sect. ^4| that the final structure de- 
pends strongly on the pattern speed Ob and a little on 
the bar strength eo. This is also evident in Fig. ||; the 
global profiles in the same Ob are quite similar, but P a in- 
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Fig. 8. The probability that an observed rotation velocity would exceed a C rit times that inferred from the galaxy 
potential. Probability P a is averaged over the time t = 300 — 500 Myr when the systems are in quasi-steady states. 
Models with different pattern speeds fib and bar strengths eo are arranged vertically and horizontally, respectively. 



creases slightly with increasing eo. Model A, B and C have 
no ILRs, not showing non-axisymmetric structures in the 
central regions (Sect. 2.4), thus the Pa-profiles are similar 
to the step functions. Model D, E and F show the most 
prominent streaming motions in their central regions, and 
therefore they have the largest P a . Although models G, 
H and I have the same fib, the P Q -profile for model G is 
different from those for models H and I, because gaseous 
X2-orbits remain in models H and I, but not in model G 



(see Sect. |2J). 

Fig. ^| shows the probability P 7 vs 7 cr ;t . All the plots 
show properties similar to those in Fig. |8[ In model E for 
example, P 7 = 0.15 at 7 cr i t = 5.0 means that the central 
galactic mass derived from an observed rotation curve is 



overestimated by a factor of five in the probability ofl5%. 
P 7 s in all models become zero at 7 cr jt = 6.0, meaning that 
the central mass from an observed rotation curve can be 
overestimated by at most a factor of six in our models. 

4. Discussion and Summary 

Based on gas dynamical calculations in a fixed galactic 
potential with a weak bar-like distortion, we estimated 
errors in mass estimation from the rotation curves, and 
calculated the probability that observations would suffer 
from such errors. We found that, as well as the final mor- 
phologies of gas disks, the probability strongly depends 
on the pattern speed of a bar f2b, and weakly on the 
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Fig. 9. The probability that an observed mass would exceed 7„it times that inferred from the galaxy potential by 
chance. Probability P 1 is averaged in t = 300 — 500 Myr when the systems are in quasi-steady states. Models with 
different pattern speeds f2b and bar strengths eg are arranged vertically and horizontally, respectively. 



bar strength e$. Among our nine models, the probabil- 
ity for the errors becomes maximal for the models with 
f2b = 0.8 X (fi — «/2) max ; the 15% of them have an error 
of a factor of five in mass estimation, if we observe the 
disks from an arbitrary viewing angle, and if we define 
rotation curves as the highest- velocity envelope of the p-v 
diagrams. Even in those erroneous cases, the galactic mass 
is not overestimated by more than a factor of six. In all of 
our models, the overestimation in mass is more probable 
than the underestimation. We consider only some partic- 
ular cases for a weak bar, thus cannot obtain general con- 
clusio ns. The above estimation however must be a guide- 



line to consider the central galactic mass derived from an 
observed rotation curve. 



Conventionally, rotation curves have been often de- 
fined as the peak-intensity velocity or intensity-weighted 
mean velocity of p-v diagrams. However, Bofuc (1996) 
pointed out that these methods underestimate the rota- 
tion velocity, particularly in the central region, because 
the finite beam size causes the confusion with the gas with 
lower line-of-sight velocities on the p-v diagram; this ef- 

~ (pOOl in their 



Koda et al 



feet is also demonstrated in 
Fig. 15). For rotation curves in highly inclined galaxies, 
this confusion can not be avoidable. Alternatively, the 
envelope-velocity of the p-v diagram is better suited to 
trace the cen tral mass distribution (sofue, 1996; Bofue & 
2001). Therefore we defined the highest-envelope 



Rubi: 



velocity as our rotation curves in the above study. Here 
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Fig. 10. Same as Fig. but rotation curves are derived by taking the density-weighted mean velocity rather than the 
most rapidly rotating envelope of the p-v diagram. 



we repeat the same analysis for a comparison, using rota- 
tion curves derived from the density- weighted mean veloc- 
ity, and shows the results in Fig. [[(j. P 7 s are always less 
than those in Fig. || and are almost zero at 7 cr i t = 1.0. 
This means that the mass derived from the mean-velocity 
rotation curves are almost always underestimated in the 
central regions of galaxies. These results suggest that the 
conventional method for deriving rotation curves from p- 
v diagrams is not also relevant to estimate the mass in 
galaxies with bar-like distortions. 



Sofuc ct al. (1999) showed that most of the rotation 
curves rise steeply from the centers, reaching high ve- 
locities of about 100 — 300 km s -1 in the innermost re- 
gions. Owing to the large fraction of the rotation curves 
with these high central velocities, they discussed the idea 



that these velocities should be attributed to massive cores 
rather than to bars. We may have a chance to statistically 
clarify whether or not the massive cores exist by compar- 
ing a probability such as ours with the observed fraction 
of rotation curves with high central velocities. When we 
define the probability P averaged in all types of barred 
and non-barred galaxies by 



P = / P a (n h ,e )f(Sl h ,e )dn h de , 



(9) 



where f(£lh, eo) is a distribution function of galaxies with 
a pattern speed f2b and bar strength eo> the existence of 
massive cores is confirmed if the fraction of galaxies with 
the central high velocities is more than P. For example, 
using our maximum calculated probability P Q | ma x = 0.4 
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and an observed fraction of barred galaxies /bar ~ 0.6 Wada, K., Sakamoto, K., & Minezaki, T. 1998, ApJ, 494, 
( Knapcn ct al. , [2000 ), P could be very roughly calculated 236 



as P < P a \ max x / bar = 0.24. Of course, we need more Wada, K., Taniguchi, Y., Habe, A. et al. 1994, ApJ, 437, 
intensive studies for a number of barred potentials and 123 
parameters, and more precise knowledge of the distribu- 
tion function of parameters for bars. 
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